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Abstract We analyze the occurrence frequency distributions of peak fluxes P, 
total fluxes E, and durations T of solar flares over the last three solar cycles 
(during 1980-2010) from hard X-ray data of HXRBS/SMM, BATSE/CGRO, 
and RHESSI. From the synthesized data we find powerlaw slopes with mean 
values of — 1.72 ± 0.08 for the peak flux, aE = 1-60 ± 0.14 for the total flux, 
and ar — 1.98±0.35 for flare durations. We find a systematic anti-correlation of 
the powerlaw slope of peak fluxes as a function of the solar cycle, varying with an 
approximate sinusoidal variation ap{t) = ao + Aacos [27r(t — io)/?cj/c;e], with a 
mean of ao = 1.73, a variation of Aa — 0.14, a solar cycle period Tcyde = 12.6 
yrs, and a cycle minimum time — 1984.1. The powerlaw slope is flattest during 
the maximum of a solar cycle, which indicates a higher magnetic complexity of 
the solar corona that leads to an overproportional rate of powerful flares. 
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1. Introduction 

In the observational part of this study (Paper I) we focus on the statistics of 
solar-flare hard X-ray fluxes during the course of the last three solar cycles, while 
theoretical modeling is referred to Paper II. Solar flares are catastrophic events 
in the solar corona, most likely caused by a magnetic instability that triggers 
a magnetic reconncction process, producing emission in almost all wavelengths. 
Since the emission mechanisms are all different in each wavelength, such as 
nonthermal bremsstrahlung (in hard X-rays), thermal bremsstrahlung and free- 
bound or recombination radiation (in soft X-rays and EUV), gyrosynchrotron 
emission (in microwaves), plasma emission (in metric and decimetric waves), 
etc. , estimates of the energy contained in each flare event strongly depends on 
the emission mechanism, and thus on the wavelength. In this study we concen- 
trate on the hard X-ray wavelength. Hard X-ray emission in solar flares mostly 
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results from thick-target bremsstrahlung of nonthermal particles accelerated in 
the corona that precipitate into the dense chromosphere. Thus, the hard X- 
ray flux is the most direct measure of the energy release rate, and thus is 
expected to characterize the energy of flare events in a most uncontaminated 
way, while emission in other wavelengths exhibit a more convolved evolution 
of secondary emission processes. The main observables that are available for 
flare statistics in hard X-rays are the peak flux P, the total flux or fluence E 
(defined as the time-integrated flux over the entire event), and the total time 
duration T of the event. In this study, we present a comprehensive compilation of 
occurrence frequency distributions of these observables obtained in hard X-rays, 
and investigate whether their behavior is different during solar cycle minima, 
including the current anomalous solar minimum. 

2. Statistics of Solar Flares in Hard X-rays 

In this Section we describe observed occurrence frequency distributions of solar 
flare hard X-ray parameters (in chronological order), discuss the properties of 
the datasets (hard X-ray energies, observational epochs, instruments), compile 
the results in Table 1 (powerlaw slopes of fluxes, fluences, durations, number of 
events, instruments, and references), and derive synthesized distributions that 
serve as reference values of the average solar flare activity. In a subsequent 
Section wc investigate whether deviations from these reference values can be 
found during the solar cycle. 

One of the earliest reports of a frequency distribution of solar hard X-ray flare 
fluxes was made by Datlowe et al. (1974), who published a frequency distribution 
of 123 flare events detected in the 20-30 keV energy range above a threshold of 
> 0.1 photons (cm~^ s""'^ keV""'^) with the OSO-7 spacecraft during the period 
of 10 October 1971 - 6 June 1972, finding a powerlaw slope of /3p ~ 0.8 for 
the cumulative frequency distribution. For compatibility we list only powerlaw 
slopes a of differential frequency distributions in Table 1, and use the conversion 
a = (3 + 1 when needed. We list also the number of events, which is a good 
indicator of the statistical uncertainty of the powerlaw slope fits (Figure 6). 

A sample of 25 microflares of smaller size were detected at 20 keV with a 
balloon-borne instrumentation of University of California Berkeley (UCB) during 
141 minutes of observations on 27 June 1980, yielding a powerlaw distribution 
with a slope of /3 ss 1 (Lin et al. 1984). 

Many more events were observed with the JJard X-Ray Burst Spectrome- 
ter (HXRBS) onboard the Solar Maximum Mission (SMM) spacecraft, which 
recorded 6775 flare events during the 1980-1985 period, exhibiting a powerlaw 
distribution of peak count rates with a slope of ap = 1.8, extending over four 
orders of magnitude (Dennis 1985). 

A next mission with hard X-ray detector capabilities was the Compton Gam- 
ma Ray Observatory (CGRO). Although it was designed to detect gamma-rays 
from astrophysical objects, it also detected solar flares systematically during 
the period of 1991-2000. Using the Burst And Source Transient Experiment 



SOLA: ms.tex; 8 October 2010; 3:07; p. 2 



Self-Organized ( 'ril iealily 



Table 1. Frequency distributions measured from solar flares in hard X-rays and gam- 
ma-rays. References: ^) Datlowe et al. (1974); 2) Lin et al. (1984); ^) Dennis (1985); "') 
Schwartz et al. (1992); 5) Crosby et al. (1993); ®) Biesecker et al. (1993); '') Biesecker et 
al. (1994); ^) Crosby (1996); ") Lu et al. (1993); 1°) Lee et al. (1993); ") Bromund et 
al. (1995); i^) Perez Enriquez and Miroshnichenko (1999); ^'^ Su et al. (2006); Christe 
et a/. (2008); Lin et a/. (2001); Tranquille et al. (2009). 
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(BATSE), statistics of flares with energies > 25 keV was sampled and more 
detailed powerlaw distributions of peak fluxes were reported with values of ap = 
1.61 ±0.03 (Schwartz et al. 1992), ap = 1.75 ±0.02 (Biesecker et al. 1993), and 
ap = 1.79 ± 0.04 (Biesecker et al. 1994) for BATSE. Biesecker (1994) noticed 
slight differences of the powerlaw slope during low activity {a.p = 1.71 ± 0.04) 
and high activity periods {ap = 1.68 ± 0.02), with the powerlaw slope usually 
flatter for high-activity periods. 

A systematic study of flares observed with HXRBS over the entire mission 
duration of 1980-1989 was conducted by Crosby et al. (1993), measuring peak 
count rates Pds (cts s~^), converted into photon fluxes Pph (photons cm~^ s~^ 
keV"""^) at energies > 25 keV, peak HXR spectrum-integrated fluxes Px (photons 
cm~i s~^), peak electron fluxes Pe (ergs s~^), flare durations T, and time- 
integrated total energies in electrons Ee (ergs), for four different time intervals 
of the solar cycle. In Table 1 we list the values for the time ranges of 1980-1982 



SOLA: ms.tex; 8 October 2010; 3:07; p. 3 



M. AsclnvaiKlcn 



(7045 events; solar maximum phase), 1983-1984 (1008 events), 1985-1987 (545 
events; solar minimum phase), and 1988-1989 (3874 events). The values of the 
powerlaw slopes change by $ 2% during different periods of the solar cycle. 

From the Wide Angle Telescope for Cosmic Hard X-Rays (WATCH) onboard 
the Russian satellite GRANAT, a sample of 1546 flare events was observed at 
energies of 10-30 keV during 1990-1992, yielding similar powerlaw slopes for 
peak count rates, ap = 1.59 ± 0.02, and flare durations ax = 2.28 ± 0.08 as 
reported before (Crosby, 1996; Crosby et al. 1998). However, it was noted that 
the frequency distribution of flare durations exhibits a gradual rollover for short 
flare durations, approaching a slope of ~ 1, so it cannot be fitted with a 
single powerlaw distribution over the entire; range of flare durations. From the 
PHEBUS instrument on GRANAT, which is sensitive to gamma-ray energies, 
Perez Enriquez and Miroshnichenko (1999) analyzed 110 high-energy solar flares 
observed in the energy range of 100 keV - 100 McV and found the following 
powerlaw slopes: ap = 1.80 ± 0.01 for (bremsstrahlung) hard X-ray fluxes at 
> 100 keV, ap = 1.38 ± 0.01 for photon energies at 0.075 - 124 MeV, ap = 
1.39 ± 0.01 for bremsstrahlung at 300 - 850 keV, as = 1.50 ± 0.03 for the 
511 keV annihilation line fluence, as = 1-39 ± 0.02 for the 2.223 MeV neutron 
line fluence, and as = 1.31 ± 0.01 for the 1-10 MeV gamma-ray line fluence. 
The flatter powerlaw slopes of the occurrence frequency distributions at higher 
energies could possibly be explained by flatter spectra (Section 3.1). 

Using data from a > 25 keV hard X-ray detector onboard the ISEE-3/ICE 
spacecraft during 24 August 1978 and 11 July 1986, Lu et al. (1993) determined 
the frequency distributions of the peak luminosity P (erg s~^), the energy E 
(erg), and flare duration T (s) and found that the measured distributions could 
be best fitted with a cellular automaton model that produced powerlaw slopes 
of ap = 1.86, as = 1.51, and ar = 1.88. The fits of the distributions included 
an exponential rollover at the upper end, which explains why they inferred a less 
steep slope for durations than previously reported. Lee et al. (1993) analyzed the 
same data and determined the correlations and frequency distribution powerlaw 
slopes with special care of truncation biases and obtained similar values for 
ISEE-3 {ap = 1.75, = 1.62, ar = 2.73) as Crosby et al. (1993) for HXRBS. 
A third study was done with the same data (Bromund et al. 1995), where the 
energy spectrum was also calculated to determine different energy parameters, 
similar to the study of Crosby et al. (1993), finding the following powerlaw 
slopes: ap = 1.86, 2.00 for the peak photon fiux Pph (photons cm~^ s~^), 
ap = 1.92, ...,2.07 for the peak electron power Pg (erg s~^), as = 1.67, 1.74 
for the total electron energy Ee (erg), and ax = 2.40, 2.94 for the total 
duration T (s), where the range of powerlaw slopes results from the choice of the 
fitting range. The flare duration T was defined at a level of 1/e times the peak 
count rate. 

From the latest solar mission with hard X-ray capabilities, the Ramaty High- 
Energy Solar Spectroscopic Imager (RHESSI) spacecraft, frequency distributions 
were determined in the 12 - 25 keV energy band from 2002 2005 (Su et al. 2006), 
finding powerlaw slopes of ap = 1.80 ± 0.02 for the peak fiuxes, and a broken 
powerlaw aT=0.9 - 3.6 for the flare duration, similar to previous flndings (e.g., 
Crosby et al. 1998). Christe et al. (2008) conducted a search of microflares and 
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Peak counts P(cts/s) 

Figure 1. Occurrence frequency distributions of hard X-ray peak count rates P (cts/s) ob- 
served with HXRBS/SMM (1980 - 1989), BATSE (1991 - 1993), and RHESSI (2002 - 2010), 
with power law fits. An average pre-flare background of 40 cts s~^ was subtracted from the 
HXRBS count rates. Note that BATSE/CGRO has larger detector areas, and thus records 
higher count rates. 

identified a total of « 25, 000 events observed with RHESSI during 2002 - 2007 
and investigated the frequency distributions at lower energies, finding powerlaw 
slopes of ap = 1.50 ± 0.03 for 3 - 6 keV peak count rates Pph (cts s~^), ap = 
1.51 ± 0.03 for 6-12 keV peak count rates, and ap = 1.58 ± 0.02 for 12 - 
25 keV peak count rates. Converting the peak count rates P into total energy 
fluxes by integrating their energy spectra, Christe et al. (2008) find an energy 
distribution with a powerlaw slope of as — 1-7 ± 0.1, with an average energy 
deposition rate of < 10^^ erg s~^. It is interesting that this microflare statistics 
is fairly consistent with overall flare statistics, even though it represents only a 
subset in the lowest energy range. 

Flare statistics was also gathered from the Solar X-ray/ Cosmic Gamma-Ray 
Burst Experiment (GRB) onboard the LOjsses spacecraft (Tranquille et al. 2009), 
finding similar results for > 25 keV events, i.e. , a powerlaw slope of ap = 
1.60 ± 0.04, which steepens to ap = 1.75 ± 0.08 if the largest events with pulse 
pile-up are excluded. 
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Figure 2. Occurrence frequency distributions of hard X-ray total counts or fiuence E (cts) 
observed with HXRBS/SMM (1980 - 1989), BATSE (1991 - 1993), and RHESSI (2002 - 2010), 
with powerlaw fits. An average pre-flare background of 40 cts s~^ muhtphed with the flare 
duration was subtracted in the total counts of HXRBS. 

We re-compiled statistics from existing flare catalogs from the three instru- 
ments HXRBS/SMM (1980 - 1989), BATSE/CGRO (1991 - 1993), and RHESSI 
(2002 - 2010) and show summary plots of the resulting frequency distribu- 
tions in Figures 1-3. The powerlaw slopes are obtained from weighted linear 
regression fits, using the Poisson statistics of the number of events in each 
logarithmic bin (see Section 3.5). The average powerlaw slope for peak fluxes 
P is ap = 1.72 ± 0.08 (Figure 1). The corresponding frequency distributions of 
total counts or fluences are shown in Figure 2, which have an average powerlaw 
slope of as — 1.60 ± 0.14 in the range oi E > 10^ cts. The distributions of flare 
durations are shown in Figure 3, which exhibit an average of ax = 1-98 ± 0.35, 
with a tendency of a rollover at the low end. Thus, our synthesized reference 
values are, 

N{P)ocP-"^ ap = 1.72 ±0.08 

N{E)o^E-'^^ = 1.60 ±0.14 (1) 

N{T)(xT-°'^ aT = 1.98 ±0.35 

which are compatible with most of the published values listed in Table 1. Thus, 
we can consider these synthesized values as representative means, averaged from 
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Total duration T(s) 

Figure 3. Occurrence frequency distributions of hard X-ray flare durations T(s) observed 
with HXRBS/SMM (1980 - 1989), BATSE (1991 - 1993), and RHESSI (2002 - 2010), with 
powerlaw fits. The flare durations for RHESSI were estimated from the time difference between 
the start and peak time, because RHESSI flare durations were determined at a lower energy 
of 12 keV (compared with 25 keV for HXRBS and BATSE), where thermal emission prolonges 
the nonthermal flare duration. 



three major missions over the last 30 years and three solar cycles, which can 
serve as a reference for the overall flare productivity of the Sun in hard X-ray 
wavelengths. 

In Figure 4 we show the correlation plots between the parameters (sampled 
in Figures 1-3) and determine linear regression fits, which yield scaling laws of 
E cx pi-36±o.04 g^^j J, ^ p0.94±o.05^ Theoretically, these correlation coefficients 
can also be calculated from the powerlaw slopes of the frequency distributions. If 
two parameters x and y have powerlaw distributions N(x) oc and N(y) cx 

y""*, and the parameters are correlated by a powerlaw relationship y (x , 
the powerlaw indices are related by /3 = {a^ — i)/{ay — 1) (e.g., Aschwanden 
2010a; Section 7.1.6). Therefore, based on the powerlaw slopes ap, as, and ax 
measured in Figures 1-3, we expect the following correlations between these 
three parameters, 

EocP^^ /3£ = (ap-l)/(ai;-l) = (1.72-l)/(1.60-l) = 1.2±0.2 
T(xP^^ = (ap - l)/(aT - 1) = (1.72 - 1)/(1.98- 1) = 0.8 ±0.2 ^ ' 
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Figure 4. Scatterplots between the total counts E(P) (left panels) or flare duration T{P) 
(right panels) versus the peak count rate P for solar flares with HXRBS/SMM (1980 - 1989) 
(top), BATSE/CGRO (1991 - 1993) (middle), and RHESSI (2002 - 2010) (bottom). Linear 
regression fits are applied for y{x) and x{y) (grey lines) and the listed values correspond to 
the average of the two linear regression fits. 
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The correlation plots shown in Figure 8, where the linear regression fit is per- 
formed for y{x) with x as independent variable, as well as for x{y) with y 
as independent variable, respectively. The averaged powerlaw slopes (of both 
linear regression fits and all three datasets) yield values of oc pi-36±o.04 ^^^^ 
T oc pO-9't±o.05^ which are fully consistent with the values derived from the 
powerlaw slopes of the frequency distributions (Eq. 2), i.e., E oc pi-2±0-2 ^^^^j 



3. Uncertainties and Errors of Powerlaw Slopes 

The determination of powerlaw slopes and their uncertainties are subject to 
methodical and data selection effects. Given the variation of observed powerlaw 
slopes as listed in Table 1, it might be useful to briefiy discuss some important 
selection effects and statistical uncertainties. 

3.1. Energy or Wavelength Bias 

In Table 1 we compare mostly fiare data in the hard X-ray regime with energies 
of > 20 keV or > 25 kcV, which all have similar powerlaw slopes of ap « 
1.6 — 1.8. Even fiares detected at lower energies > 8 keV (Lin et al. 2001) or 

> 12 keV (Su et al. 2006; Christe et al. 2008) have similar powerlaw slopes. On 
the other end, flares detected at high energies of > 100 kcV (Perez Enriqucz 
and Miroshnichenko 1999) have also similar powerlaw slopes. However, the same 
subset of flares observed at 0.075 - 124 MeV, from bremsstrahlung at 300 - 850 
keV, from the 511 keV annihilation line flucncc, from the 2.223 McV neutron 
line fluence, and from the 1-10 MeV gamma-ray line fluence, have significantly 
flatter powerlaw slopes {ap w 1.3 — 1.4). This can be explained by the flatter 
energy spectra of large flares, which causes relatively higher fluxes at higher 
energies, and thus flatter occurrence frequency distributions. Dennis (1985) has 
shown a systematic tendency of the hard X-ray spectrum as a function of the 
hard X-ray count rate, being generally flatter for large flares with high count 
rates, which predicts also flatter powerlaw slopes of the occurrence frequency 
distributions at higher energies. Similar wavelength bias effects occur also for 
flare fluxes detected in soft X-rays and EUV wavelengths, compared with hard 
X-ray fluxes, which depend on the physical emission mechanism. Therefore we 
confine our comparisons in Table 1 only to hard X-ray energies in the range of 

> 8 keV to > 100 keV here. 

3.2. Effects of Data Gaps and Spacecraft Orbits 

Data gaps due to spacecraft night or South Atlantic Anomaly (SAA) passages 
can also cause systematic biases, because they tend to shorten fiare durations 
(due to interruption) and underestiamte the total flare counts (due to incomplete 
sampling). These effects can in principle be simulated and their systematic effects 
on frequency distributions can be quantified this way. Numerical simulations of 
such selection and instrumental effects are also discussed in the context of waiting 
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Peak counts P(cts/s) 

Figure 5. Occurrence frequency distributions of hard X-ray peak count rates P (cts/s) ob- 
served with HXRBS/SMM (1980 - 1989), without bacl<ground subtraction (ap = 2.06 ±0.05; 
blue), and with subtraction of an average background of B = 40 cts s"! {ap = 1.73 ± 0.01; 
red). 

time distributions {e.g. , Aschwanden and McTiernan 2010). The distribution 
of peak fluxes is less affected by data gaps, as long as the peak of a flare is 
observed. Visual screening of flare light curves was carried out for events detected 
with HXRBS/SMM, so that missed flare peaks could be identified in the flare 
catalog (Dennis, private communication). For RHESSI flare catalogs, however, 
no correction for spacecraft night datagaps has been applied to flare peak fluxes 
and fluences so far, and long-duration flares are counted multiple times for each 
spacecraft orbit (McTiernan, private communication). 

3.3. Pre-Event Background Subtraction 

Most flare catalogs list the peak counts of flares without preflare background 
subtraction, which has little impact on large flares, but leads to a strong over- 
estimate of fluxes for the samllest flares, and thus systematically overestimates 
the powerlaw slope. We demonstrate this for the dataset of HXRBS/SMM flares. 
Figure 5 shows two (weighted) powerlaw flts, a steeper slope of ap = 2.06 ±0.05 
for the uncorrected flare peak counts, and a flatter slope of ap — 1.73±0.04 after 
subtraction of an estimated average background of 40 cts s~^. This correction 



SOLA: ms.tex; 8 October 2010; 3:07; p. 10 



Self-Organized ( 'ril iealily 



was also applied in the original analysis in Crosby et al. (1993) and we use it here 
too. Note that the background-uncorrected distribution shows also a steepening 
of the powcrlaw slope at the lower end. 

For the RHESSI flare catalog, the preflare background is estimated to be 
« 10 — 20 cts s~^ per detector, the low limit applying to low spacecraft latitude, 
the high limit for high latitutes, with a mean of » 12 cts per detector (Jim 
McTiennan, private communication). 

The subtraction of an average preflare background value is only a first-order 
correction, which may cause even events with negative peak fluxes when the 
background is time-varying. The optimum method would be to measure the 
preflare background for each event separately and to subtract it from the peak 
fluxes, as well as from the fluences before time integration. 

3.4. Instrument Sensitivity 

The instrumental sensitivity usually translates into a fixed flux threshold, which 
causes a sharp lower cutoff for peak flux distributions (Figures 1 and 5), but a 
gradual roll-over for total fluence and duration of events (Figure 2 and 3) . The 
roll-over limits the scale-free range of the distribution over which a powerlaw 
can be fitted, and moreover produces a systematic bias for flatter powerlaw 
slopes when part of the roll-over is included in the powcrlaw fit. If powerlaw 
fits are weighted by the number of events per bin in a linear regression fit to 
a histogram, the lowest bins will have the hightest weight, and thus a gradual 
roll-over at the lower end will produce a too fiat powerlaw slope, when compared 
with the fit in the upper bins of the distribution. We determine the lower bound 
of the powerlaw fit range by the x^-criterion, which reliably detects a roll-over 
by an excessive deviation from a powerlaw fit in terms of the expected standard 
deviation. 

3.5. Linear Regression Fit Methods 

A common method to determine powerlaw slopes is a linear regression fit to 
a histogram in a log{N) — log{S) diagram, which can be done in two ways. If 
the frequency distribution closely follows a powerlaw function, a weighted linear 
regression fit can be performed in log-log space, using the number of events per 
bin as relative weights. However, if the frequency distribution exhibits significant 
deviations from a straight powcrlaw (e.g. , Figures 2 and 3), a weighted powerlaw 
fit can only be performed in a powerlaw-like subinterval. Often an unweighted fit 
(with equal weight in log-log space) is carried out in such cases in order to obtain 
an approximate powerlaw slope. For the results given in literature (as compiled 
in Table 1), it is not always clear whether the authors performed a weighted or 
unweighted linear regression fit, which may significantly deviate from each other 
depending on the fitted interval. In this study here, we perform only weighted 
fits in well-defined intervals as described below. 

In a log-log histogram, the bins Aa;^ are chosed to be equidistant on a log- 
scale, but they are exponentially increasing on a linear scale. If N{x) w x~°' is 
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the fitted powerlaw function, the number of events Ni per bin is, 

fXi+Axi 



J Xi 



N{x)dx w N{xi)Axi . 



(3) 



The statistical uncertainty a|°^ in bin Axi due to Poisson statistics is then on a 
log-scale, 

= log(iVi + ^/Nl) - log(A^i) , 



(4) 

and the statistical weighting factor Wi for a linear regression fit that minimizes 
the squared standard deviations in log-log space is then. 



Wi 



■ 



(5) 



The linear regression fit in log-log space finds then the following best-fit powerlaw 
function, 



logiNp') = No + alog{x,) 



(6) 



A criterion for the goodness of a fit is the x^-tcst, expressed with the so-called 
reduced x^-va_/ue, which expresses the mean deviation of the data points from 
the fitted function in units of expected standard deviations. 



Xred 



1 ^{Ni-Nfy 



-,1/2 



(n - Ufree) ^ (af »)2 
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where n is the number of datapoints (or histogram bins here) and n/^ee is the 
number of free parameters of the fitted function, which is n free = 2 for a linear 
regression fit (with free parameters A^o a-nd a in Eq. 6). However, since the upper 
part of the histogram in log-log distributions generally include bins with zero or 
very few events, the validity of Poisson statistics v^/Vj" is hampered for Ni < 10 
and has to be replaced by C-statistics (Cash 1979), which has the following 
expression corresponding to the reduced x^-criterion. 



Xcash — 



(n - Ufree) 



Niln{Ni/Np')-{Ni-NP') 



i=l 



1/2 



(8) 



where A"/** is the theoretical expectation value based on the best-fit powerlaw 
function and has to be positive for all fitted bins. For an acceptable powerlaw 
fit, a value of xcash ~ 1-0 is expected. The linear regression fit yields then a 
powerlaw slope a and the probable uncertainty <Ja due to Poisson statistics, or 
C-statistics, respectively. 

3.6. Monte-Carlo Simulations of Powerlaw Uncertainty 

An alternative method to infer the uncertainty of powerlaw slopes is a Monte- 
Carlo simulation of generating a powerlaw distribution. Mathematically, a de- 
sired distribution function f{x) can be generated from uniformly distributed 
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Figure 6. Monte-Carlo simulations of powerlaw distributions are carried out for a theoretical 
powerlaw slope of a = 1.75 with different random representations a<;cording to Equation 11. 
A single distribution based on AT = 1000 events is shown with powerlaw fit a = 1.73 (top 
left panel). The Gaussian distribution of powerlaw slopes obtained from Nrun = 1000 runs 
with different random representations yields a mean of a = 1.74 ± 0.03 (top right). These 
Monte-Carlo siiuulatioiis are repeated for various iminber of events N = 10^, .... 10''', which 
show a convergence to the mean value and standard deviations towards the theoretical value 
a = 1.75 (bottom left). The scaling of the uncertaingy Sa{N) as a function of the number of 
events N is also shown (bottom right). 



random values p in the range of [0,1] by using the inverse function F~^{x) 
of the integral function F{x) of f{x). Here we wish to generate a frequency 
distribution of energies E in the form of a powerlaw function p{E) with slope a 
(Aschwanden 2010a, Section 7.1.5), 



p{E) = {a- 1)E- 



(9) 



which fulfills the normalization p{E)dE = 1. The total probability p{E) in 
the range [0, E] is then the integral function of p{E) (Eq. 9), 



p{E) = / p{e)de = (a - l)e-"de 
Jo Jo 



1 _ _e;(i-") 



(10) 
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The inverse function E{p) of p{E) (Eq. 10) is 

£;(p) = [i-p]^/(^-") . (11) 

In Figure 6 (top left) we use a random generator that produces 1000 values pi 
uniformly distributed in the range of [0, 1], choose a powerlaw index of a = 1.75, 
and use the transform Eq. (11) to generate values Ei = [1 — pj]^/^^"") and sample 
the frequency distribution of the 1000 values Ei, which produces a powerlaw 
function p{E) = [a — 1)E~°' as defined in Eq. (9), with a best-fit slope of 
a = 1.73. 

We repeat the generation of powerlaw histograms with = 1000 different 
random representations Pi and find for the powerlaw slopes a Gaussian dis- 
tribution with a mean and standard deviation of a = 1.74 ± 0.03 (Figure 6, 
top right). The spread of a- values is expected to depend on the number of 
events in a distribution. In order to quantify this dependence we repeat the same 
exercise for different numbers of events in the range of N = 10^, 10^, which is 
shown in Figure 6 (bottom left). The standard deviations of the powerlaw slope 
distributions is found to depend on the number of events (per distribution) as 
da « 0.45/log(A'')^ (Figure 6, bottom right). Repeating the same Monte-Carlo 
simulations for a range of powerlaw slopes (a = 1.5, 2.0) we find the following 
general dependence, 

6aia,N)^ia-0.75)^. (12) 

This uncertainty of the powerlaw slope obtained from Monte-Carlo simulations 
is typically a factor of « 1.5 ± 0.4 larger than the formal uncertainty obtained 
from linear regression fits (according to the statistics given in Table 2). Both 
methods are based on the Poisson statistics of random numbers and are expected 
to lead to commensurate uncertainties of powerlaw slopes. However, while the 
uncertainty in a linear regression fit is based on a single dataset, the uncertainty 
from Monte- Carlo simulations includes in addition the scatter of many random 
representations and thus yields a slightly higher uncertainty. Thus, we consider 
the latter method as a more reliable measure or uncertainties when comparing 
the best-fit powerlaw slopes among two different datasets. 



4. Solar Flare Statistics versus Solar Cycle 

We investigate now whether the statistical distributions of the most reliably 
determined flare parameter, namely the peak flux P, exhibits significant changes 
during the solar cycles. From the time profiles of the monthly flare rates shown 
for solar cycles 21, 22, and 23 (Figure 5), we identify approximately the following 
time intervals for the solar minima: 1985 - 1987, 1995 - 1997, and 2007 - 2010. 

For the first solar minimum we have coverage by HXRBS/SMM data. We 
break down the dataset into 2-year intervals and show the frequency distributions 
of the peak flux separately for 1980-1981, 1982-1983, 1984-1985, 1986-1987, and 
1988-1989 in Figure 8 (left column). We perform weighted linear regression fits in 
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Figure 7. Top: Monthly averages of solar flare rates observed during the last three solar 
cycles in hard X-rays with HXRBS/SMM (1980 - 1989), BATSE/CGRO (1991 - 2000), and 
RHESSI (2002 — 2010). Bottom: Monthly averages of solar flare rate observed in soft X-rays 
with GOES, including events above the C3-class level. 

the lower half of the distributions, bracketed by the values Pi < P < P2 (dotted 
lines in Figure 8), where Pi corresponds to the next bin above the maximum of 
the frequency distribution, and P2 is chosen to be higher by a constant factor of 
20. The number of events n, the lower and upper bounds Pi and P2, the powerlaw 
fits ap with the formal fit uncertainty aa, the uncertainty 6ap obtained from 
the Monte-Carlo simulations (Equation 12), and the reduced x^-value according 
to C-statistics (xcash) sue all listed in Table 2 and shown in Figure 8. We arc 
doing the same procedure for the BATSE/CGRO dataset by breaking it down 
into the two intervals of 1991 and 1992-1993, which covers part of Solar Cycle 22. 
In addition, we perform the same procedure for RHESSI data by breaking down 
the dataset into the following 2-year intervals: 2002-2003, 2004-2005, 2006-2007, 
2008-2009, and a 1-year interval of 2010. Thus, we determined the powerlaw 
slope in the lower half of the frequency distribution for a number of 12 time 
intervals with a length of 1-2 years during the three different solar cycles. 

In order to investigate the time dependence of the powerlaw slopes of peak 
flux distributions we plot the values ap±Sap in the middle panel of Figure 9, 
which exhibits a sinusoidal variation diiring the solar cycles. For comparison, we 
plot also the published values of datasets with a timespan less than a half solar 
cycle (< 6 years) and a minimum of n > 100 events, which includes datasets 1, 
3-8, and 12-14 from Table 1 (Figure 9, top panel). The datasets labeled with 
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Figure 8. Occurrence frequency distributions of hard X-ray peak count rates P (cts/s) ob- 
served witii HXRBS/SMM, BATSE/CGRO, and RHESSI, broken down into 2-ycar intervals. 
The powrerlaw fits virere determined with weighted linear regression fits applied to the lower 
half of the sampled distributions , within a range ratio of P2/P1 = 20 (dotted vertical lines 
and thick lincstylc for best-fit power law functions). The best-fit values are also listed in Table 
2 and plotted as a function of time in Figure 9. 
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Table 2. Frequency distributions measured from solar flares in hard X-rays during the solar 
cycle. 



Instrument 


Years 


Number 


Lower 


Upper 


Powerlaw 


Goodness 






of 


fit 


fit 


slope of 


of fit 






events 


bound 


bound 


peak counts 








n 


Pi 


P2 


Qp (Sap) 


\,Cash 


HXRBS 


1980-1981 


4460 


43 


848 


1.77±0.02 (±0.05) 


1.43 


HXRBS 


1982-1983 


2642 


43 


848 


1.77±0.03 (±0.05) 


0.66 


HXRBS 


1984-1985 


581 


43 


818 


1.91±0.06 (±0.08) 


1.16 


HXRBS 


1986-1987 


437 


43 


8-18 


1.75±0.08 (±0.07) 


1.03 


HXRBS 


1988-1989 


3286 


31 


439 


1.60±0.03 (±0.04) 


1.34 


BATSE 


1991-1991 


2389 


610 


11787 


1.61±0.03 (±0.05) 


0.56 


BATSE 


1992-1993 


1720 


610 


11787 


1.69±0.03 (±0.05) 


0.90 


RHESSI 


2002-2003 


6678 


31 


439 


1.59±0.02 (±0.04) 


1.69 


RHESSI 


2004-2005 


3930 


31 


439 


1.66±0.03 (±0.04) 


1.78 


RHESSI 


2006-2007 


958 


31 


439 


1.83±0.06 (±0.07) 


1.31 


RHESSI 


2008-2009 


581 


31 


610 


1.83±0.05 (±0.07] 


0.65 


RHESSI 


2010-2010 


1496 


43 


848 


1.89±0.04 (±0.06] 


0.79 



17-19 in the middle panel correspond to HXRBS, BATSE, and RHESSI as listed 
in Table 2 and analyzed in Figures 1-4 and 8. 

We are now characterizing the time variation of the powerlaw slope ap{t) as 
shown in Figure 9. First we test the hypothesis wliether the data are consis- 
tent with a constant value, using simply the x^-criterion (Equation 7) for the 
expectation value of a constant mean value (horizontal lines in Figure 9). For 
the published values we find a mean of ap = 1.72 ± 0.08 and a Xred = 2.01 
(Figure 9 top panel), and similarly for our re- analyzed powerlaw slopes with 
weighted linear regression fits in the lower half of the frequency distributions, 
i.e. ap = 1.74 ± 0.11 and Xred = 2.7 (Figure 9 middle panel). Based on this 
unacceptable goodness-of-fit value Xred ^ 2 we can reject the hypothesis that 
the powerlaw index is constant during the solar cycles. 

Next we test the hypothesis whether the powerlaw slope ap{t) varies with a 
sinusoidal variation, 

ap{t) =ao + Aa cos [2n{t - toj/T^yde] ■ (13) 

For the published values, using the uncertainties that are quoted in the publi- 
cations, we find no acceptable fit, based on a goodness-of-fit value of Xred = 2.5 
(Figure 9, top panel). The discrepancy comes mostly from RHESSI data, where 
two quite different values were measured during the same time interval, namely 
ap = 1.80 ± 0.02 (dataset 13; Su et al. 2006) versus ap = 1.58 ± 0.02 (dataset 
14; Christe et al. 2008). 

For the powerlaw fits obtained in this study, however, we find an acceptable 
fit with Xred = 1-15 and best-fit parameters of: mean powerlaw slope ap = 1.73, 
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Figure 9. The powerlaw slope ap of the peak flux P as a function of time during the last 
three solar cycles. The diamond symbols and numbers indicate the mean values cip and the 
corresponding reference number quoted in Table 1. The horizontal lines indicate the time 
intervals of observation and the vertical bars indicate the cited uncertanties. The functional 
dependence of the powerlaw slope ap(t) on the solar cycle is tested with a fit of a constant 
value (horizontal lines) and with a sinusoidal curve (dashed curve in top panel, and thick curve 
in middle panel). The bottom panel shows the GOES flare rate for comparison (see also Figure 
7). 



amplitude of variation Aa = 0.14, solar cycle period Tcyde = 12.6 yrs, and a 
cycle minimum time to = 1984.1 (Figure 9, middle). This is a novel result that 
shows for the first time the solar cycle dependence of the powerlaw slope of flare 
hard X-ray count rate distributions. For comparison we plot the monthly GOES 
flare rate (for > C3 class) in the bottom panel of Figure 9 (also shown in Figure 
7). The variation of the powerlaw slope is clearly anti-correlated to the flare rate 
or sunspot number, bcin steepest during solar minima and flattest during solar 
maxima. The best-fit solar cycle period is Tcyde = 12.6 years applies to the mean 
behavior during the fitted period of 1980-2010, but is heavily weighted by the 
longest last Cycle 23. The actual lengths of solar cycles published by NOAA are 
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10.3 years for Cycle 21 (1976-1986), 9.7 years for Cycle 22 (1986-1996), and 12.6 
years for Cycle 23 (1996-2006). 

5. Discussion and Conclusions 

Frequency distributions of solar-flare parameters have been reported over the 
last three solar cycles from different instruments and different phases of the 
solar cycle, covering a scattered range of values (Table 1). Power law slopes were 
reported in a range of ap « 1.6 — 2.0 for peak fluxes (counts), as ~ 1.4 — 1.7 
for total fluxes (fluences), and ar ~ 1-9 — 2.8 for total durations, if we restrict 
to the hard X-ray range of w 10 — 100 keV. Averaging the powerlaw slopes from 
the largest datasets spanning entire solar cycles (using data from HXRBS/SMM, 
BATSE/CGRO, RHESSI), wc find mean values of ap = 1.72 ± 0.08 for peak 
fluxes (Figure 1), a_E = 1.60 ±0.14 for fluences (Figure 2), and ut = 1.98 ±0.35 
for flare durations (Figure 3). Investigating variations of the flare statistics on 
time scales of 1-2 years, however, we do find evidence for a significant variation 
of the powerlaw slope of the frequency distributions during the solar cycle. We 
flnd that the the powerlaw slope varies sinusoidally with an amplitude of ±8% 
during the solar cycle, in anti-correlation to the sunspot number, being flattest 
during the maximum of the solar cycle. 

Why is our obtained result of the anti-correlation of the powerlaw slope of flare 
count rates with the solar cycle not evident from previously piiblishcd values? 
First of all, the datasets have been subject to inconsistent analysis methods 
by different authors, differing in the type of linear regression flts (weighted vs. 
unweighted), the fitted range [Pv,P2\i selection effects of events (flare catalogs 
vs. automated microflare search algorithms), background treatment (with or 
without background subtraction), as discussed in Section 3. In order to identify 
some systematic effects we varied the fitting range [Pi, P2] between one to three 
decades and found that the sinusoidal variation shows up clearest for a factor 
of P2/P1 = 20, which entails the lower half of frequency distributions, typically 
GOES C-class to M-class flares, but excludes the large GOES X-class flares that 
make up the upper half of the frequency distributions. Obviously, X-class flares 
tend to appear more irregularly during the solar cycles, while the more frequently 
C-class and M-class flares vary in a better synchronized way with the solar cycle. 
The fact that the Sun produces flatter powerlaw slopes during periods of high 
flare activity indicates that the magnetic complexity is higher during activity 
periods and produces ovcr-proportionally energetic flares during these periods. 
The effect is analog to the soft-hard-soft evolution of energy spectra known 
during single flares. 

Marginal variations of the powerlaw slope with flare activity have been noted 
before. Crosby et al. (1993; Fig. 12 therein) found variations of the powerlaw 
slope during the solar cycle in excess of the statistical uncertainties of the linear 
regression fits, but no clear functional dependence on the solar cycle was foimd. 
Biesecker (1994) noticed slight differences of the powerlaw slope during low 
activity {ap = 1.71 ± 0.04) and high activity periods {ap = 1.68 ± 0.02), with 
the powerlaw slope usually flatter for high-activity periods, which is analog to 
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our finding of a flatter slope during solar cycle maxima. Bai (1993) devised a 
special maximum likelihood method to determine the powerlaw slope of flare 
distributions from HXRBS/SMM and found some variation of the powerlaw 
slope correlated with a 154-day periodicity of flare rates. He found that the size 
distrtibutions are steeper during the maximum years of Solar Cycle 21 (1980 and 
1981) than in the declining phase (1982-1984). The cycle effect is opposite to our 
findings, but it depends sensitively on the threshold or fitted range. However, 
the three studies (Bai 1993; Biesecker 1994; and this study) agree in the result 
that a flatter powerlaw slope is correlated with a higher flare activity, either on 
periods of recurrent flare activity (153.8 days) or solar cycles (12.6 years). 

What docs oiu- novel result suggest for the physical interpretation? A flat- 
ter frequency distribution implies an overproportional amount of larger events, 
compared with an averaged distribution. This means that the physical conditions 
are different during the solar cycle. One possible explanation could be that the 
magnetic complexity increases during the solar maximum, which produces more 
stressing of magnetic flelds and larger energy releases. The magnetic flelds during 
the maximum of the solar cycle are dominated by the toroidal component of 
the solar dynamo, which entails higher magnetic stresses. The magnetic field 
during the minimum of the solar cycle is simpler and is dominated by the 
dipolar poloidal field that spans from the north to the south pole. It appears that 
the relative amount of energy released in flares is modulated by the magnetic 
field complexity of the solar dynamo. Theoretical models that reproduce the 
statistical distributions of solar flares in terms of the concept of self-organized 
criticality are discussed in Paper II (Aschwanden 2010b). 
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